Perfect fluid flow from granular jet impact 
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Experiments on the impact of a densely-packed jet of non-cohesive grains onto a fixed target show 
that the impact produces an ejecta sheet comprised of particles in collimated motion. The ejecta 
sheet leaves the target at a well-defined angle whose value agrees quantitatively with the sheet angle 
produced by water jet impact. Motivated by these experiments, we examine the idealized problem 
of dense granular jet impact onto a frictionless target in two dimensions. Numerical results for the 
velocity and pressure fields within the granular jet agree quantitatively with predictions from an 
exact solution for 2D perfect-fluid impact. This correspondence demonstrates that the continuum 
limit controlling the coherent collective motion in dense granular impact is Euler flow. 



Impact is a powerful tool in physics. We use scattering 
techniques to deduce the structure of matter by analyzing 
the ejecta distribution produced by the impact of a highly 
collimated beam against a target. Here we examine a 
macroscopic version of beam impact: the collision of a 
jet of macroscopic grains against a fixed target. We find 
that to leading order, the impact of a dense stream of 
grains is described by Euler flow. That is, the granular 
dynamics simplifies greatly and is governed primarily by 
incompressibility and momentum conservation. 

Our work is motivated by experiments showing that 
the collision of a dense granular jet with a fixed target 
produces a liquid-like dynamical response pQ . At low vol- 
ume fractions, impact of a granular stream onto a fixed 
target produces a shock which quantitatively agree with 
the Mach cone produced by supersonic gas flow [2H1]. 
At intermediate volume fractions, experiments and sim- 
ulations show coexistence of a shock-like front coexist- 
ing with a dead zone [5HH]. In contrast, when the jet 
is densely-packed, the amount of scatter is significantly 
reduced. Instead, the angles of recoil are narrowly dis- 
tributed about a mean value \I/o (Fig.[l](a)). As a result, 
the granular jet is ejected away from the target in a thin 
axisymmetric sheet forming the shape of a hollow cone, 
reminiscent of the "water bell" [9] . Surprisingly, the vari- 
ation in the ejecta angle \I/o as a function of the geometric 
ratio D tar /Dj et in the granular jet system quantitatively 
tracks the variation in the water jet near the target pQ, 
where D taT is the target width and Dj et is the jet width. 

In the regime we focus on here, the volume fraction of 
grains within the jet is close to random close packing, im- 
pact creates a dead zone comprised of grains which are ef- 
fectively immobile directly over the target (see Fig.[T](b)). 
Note that because of the high packing fraction, the dead- 
zone is distinct from a shock because the packing fraction 
remains at random close packing inside and outside of the 
deadzone. 

In an earlier studies [TQl [11] , we found that removing 
frictional effects at the target eliminates the dead zone 
but leaves the ejecta dynamics essentially unaffected. Be- 
cause the existence of the dead zone is not essential to the 
emergence of the liquid-like ejecta dynamics, we here opt 




FIG. 1. (Color online) (a) Snapshot of the axisymmetric, 
hollow ejecta sheet formed by dense non-cohesive granular jet 
impact. The sheet leaves the target with an angle (from 
PQ). (b) Snapshot of a 2D simulation of granular impact onto 
a stick target: once a grain touches the target, it remains mo- 
tionless for the rest of the simulation. Here, D gra in = Dtar/25. 
A dead zone marked by the blue grains is present. The dead 
zone is distinct from a shock because the density is continuous 
across the dead zone surface. Smaller D gra i n /D ta r increases 
the coherence of the ejecta. (c) Snapshot of a 2D simula- 
tion of granular jet impact onto a frictionless, or slip, target, 
with D gra i n — L)tar/25. In 2D, the collimated ejecta form two 
outgoing jets instead of the conical shape we see in 3D. The 
collimated ejecta forms in 2D and 3D, and with or without 
the interior dead zone. 



to analyze the idealized problem of impact onto a fric- 
tionless target. The absence of a dead zone simplifies the 
dynamics, allowing us to precisely identify the continuum 
dynamics controlling the collimated ejecta. Specifically, 
we compare predictions from an exact solution for the 
impact of a perfect-fluid jet with results from a full dis- 
crete element method simulation. The two sets of results 
agree quantitatively, demonstrating that dense granular 
impact simply corresponds to Euler flow at leading order. 

In the rest of the paper, we first describe the idealized 
impact problem in more detail, then give the numerical 
results showing that the granular jet during impact re- 
mains "cold," where the energy contained in fluctuations 
is small, and that the large-scale motion produced by im- 
pact is to a high degree incompressible and irrotational. 
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These three features suggest that the impact in essence 
corresponds to the impact of a perfect fluid. Finally we 
check this idea explicitly by comparing the numerical re- 
sults of kinematic and stress fields against exact solutions 
for perfect fluid flow and obtain good agreement. Here 
we present only results from a 2D impact onto a fric- 
tionless target where £>j e t? the width of the granular jet, 
is equal to the target width D tar . Impact in 3D or for 
different D tav /D- ]et ratios do not significantly change the 
results described here. 

Previous experiments used copper and glass beads, 
both of which are fairly rigid, and found that changing 
the impact speed does not change the eject a dynamics, 
only the characteristic scale of the forces present. In our 
simulations we therefore idealize the particles within the 
granular jet as perfectly rigid circular disks. To avoid 
spurious crystallization in the densely packed jet, we 
use poly-disperse particles, with diameters uniformly dis- 
tributed between 0.8D gra i n and 1.2D gra i n . Unless other- 
wise specified, our simulations used D ]et = 200D gra i n , a 
ratio within the range used in the experiment. We have 
varied both the coefficient of friction and the coefficient 
of restitution between grains. We found, in agreement 
within the range accessed by experiments pQ that their 
effects are small. The results presented here use a co- 
efficient of friction of 0.2 and a restitution coefficient of 
0.9. 

To simulate the impact, we use an time step driven 
scheme modified to handle dense flows efficiently [T2] . 
We have compared results from the simulation against 
the experiment for 3D impact onto a frictional target and 
verified that we obtain quantitative agreement between 
the two [10]. 

In 2D, collision with the target creates two collimated 
ejecta streams (see Fig. [ijc)). The simulations with a 
frictionless target also clearly show that the entire inte- 
rior of the jet is in continual motion, with the veloc- 
ity varying smoothly from the target to the free sur- 
face. To characterize this collective motion, we time av- 
erage the motions of the individual grains over a time 
scale 30D t ar/^o spatially bin grains into bins of lin- 
ear dimension 2.5D gra i n to obtain (u(x)) and (p(x)) = 
{{? xx + a yy)/2)i velocity and pressure fields that vary 
spatially. 

We can gauge whether the collimated flow in granular 
jet impact is due to the presence of equilibrium hydro- 
dynamics by considering the energy fluctuations within 
the jet. For example, in the case of liquid jet impact un- 
der the laboratory conditions in [9], the energy density 
of fluctuations Nk^T/V normalized by the impact en- 
ergy density scale Pjet^o are about 10 3 . In this case, we 
know that the water is near thermodynamic equilibrium 
and the macroscopic flow is thought of as a small pertur- 
bation to an equilibrium liquid phase. The equilibrium 
status for granular flow is quite different from this, as we 
show here. To assess the relative importance of energy 




FIG. 2. (Color online) (a) Fluctuations in the kinetic en- 
ergy density, or the "granular temperature" Tq, divided by 
the impact energy density scale pjetUo for half of the jet near 
the target. The granular temperature Tq is about 10 ~ 3 of 
the impact energy when Dt ar / D gra i n = 200. (b) Normalized 
average granular temperature in the warm region above the 
target (+) and in the hottest region at the target edge (JC) 
plotted against the relative grain size D gTa i n / D taT . By de- 
creasing the grain size D gra i n / D ta r, we suppress fluctuations 
within the jet. 

fluctuations in dense granular jet impact, in Fig. |2](a) 
we plot the amplitude of the fluctuations in the energy 
density within the jet Tq = pj et ((u 2 } — (?/) 2 )/2, known 
as the "granular temperature," normalized by the impact 
energy density scale Pjet^o • The normalized amplitude of 
the energy fluctuations normalized by the impact energy 
scale is about 10 -3 when D tar /D gra i n = 200. In Fig.[2](b) 
we see that as the grain size D gra i n /D tar is decreased, 
the granular temperature decreases. In the continuum 
limit of infinitesimally small grains, the energy fluctua- 
tions appear to vanish. We conclude that the granular 
jet is far from equilibrium, and that the continuum flow, 
and therefore also the collimated ejecta, is not the result 
of an equilibrium hydrodynamical phase on impact. 

In Fig. [3](a) we show the non-dimensionalized diver- 
gence |V • (u}\/(Uo/D tar ) throughout the jet. We find 
that the divergence of the mean velocity field is small, 
on the order of 0.1% of the characteristic scale of the 
velocity gradient. In Fig. [3](b) we examine the non- 
dimensionalized vorticity, |V A (u)\/ (Uo/ D tar ). The nor- 
malized vorticity in the majority of the jet is, at roughly 
1% the magnitude of the velocity gradient created by the 
impact, relatively small. The only exceptions to the flow 
being incompressible and irrotational are regions near the 
free surface of the ejecta sheet, where the particles have 
ceased colliding with each other. This suggests that the 
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FIG. 3. We plot (a) the nondimensionalized divergence 
Aar(V • (u))/U and (b) vorticity Aar|V A (u)\/U for half 
of the jet. The divergence and more noticeably the vorticity 
become large on the free streams behind the target, in the 
regime of ballistic motion where grain-grain interactions are 
weak. But near the target, where the collimated response is 
produced, both the vorticity and divergence remain small. 

average motion is, to a high degree, incompressible and 
irrotational throughout the impact region. 

A perfect fluid describes the flow of an incompress- 
ible continuum material uninfluenced by thermodynam- 
ics. The vorticity of a perfect fluid does not change in 
time, hence a perfect fluid in jet impact will be irrota- 
tional before impact and remain irrotational. The finding 
that granular jet impact is largely incompressible and ir- 
rotational suggests that the continuum limit for dense 
granular impact can be idealized as a perfect fluid. 

The impact of a perfect fluid satisfies the Euler equa- 
tions 

p jet (<9 £ + u • V)u = -Vp, V • u = (1) 

everywhere in the jet interior. The surface shape of the 
jet is not known a priori, but is determined along with the 
rest of the flow when subject to the following boundary 
conditions. The pressure remains constant along the free 
surface of jet. Far upstream and downstream of the tar- 
get, the velocity within the jet is uniformally translating. 
In 2D and in steady state, an exact solution describing 
impact of a perfect fluid jet can be obtained via classical 
conformal mapping techniques [I3j [14] . 

Perfect fluid flow conserves energy, yet the granular 
jet has dissipation by means of inelasticity and Coulom- 
bic friction. For this reason, we expect and indeed find 
some discrepancy between the Euler flow solution and 
the granular jet simulation results. Specifically, because 
of these dissipative processes the velocity in the ejecta 
sheet does not recover the incident impact speed /To, as 
it should if energy is conserved. Instead the exit speed 
is reduced slightly from the incident value. Since the 




FIG. 4. (Color online) Comparison of (a) the normalized 
speed field |u(x)|/C7o and (b) the normalized pressure field 
p(x)/_p max between an exact solution for 2D perfect fluid jet 
impact (dashed lines) against a simulated 2D jet impact (solid 
lines) for half of the jet. The maximum pressure p ma x is 
achieved at the target center. Values from the simulation 
are time averaged. Here D t ar/^jet = 1- The y coordinate is 
the distance above the target along the direction of incoming 
jet (Fig. [ljc)). Euler flow captures the speed and pressure 
fields of the granular jet to leading order. 



flow is incompressible, this causes the ejecta sheet in the 
granular simulation to be slightly wider than the ejecta 
sheet created by perfect fluid impact. The same reason 
causes the ejecta angle to be slightly lower for the gran- 
ular jet impact (\I/o = 51°) than that of a perfect fluid 
(^o = 54°). 

In Fig. |4](a) we plot speed contours predicted by the 
exact solution against numerical results at 3 different con- 
tours, as well as a comparison along the center line of the 
jet. We normalize the speed by dividing the local velocity 
by the impact speed Uq. The largest discrepancy occurs 
along the boundary of the jet, as evidenced by the 0.9 
speed contour intersecting the jet surface profile from the 
numerics. This is because the velocity along the surface 
of granular jet slows slightly due to energy dissipation. 
Near the target, the contours agree quantitatively. Over- 
all the agreement is surprisingly good. 

Having found that the Euler flow solution reproduces, 
at leading order, the kinematics of dense granular impact, 
we next examine its prediction about the stress distribu- 
tion within the jet. Since Euler flow is purely driven by 
gradient of the pressure field, i.e. it supports tangen- 
tial stresses, in Fig. |4](b) we compare the pressure from 
Euler flow against the pressure from the granular impact 
simulation. We normalize the pressure by the maximum 



4 



pressure p max achieved at the center of the target. The 
pressure data from numerics agrees quantitatively with 
predictions from perfect fluid flow. However, while the 
speed contours exhibit little fluctuation throughout the 
jet interior, the pressure contours show larger and larger 
fluctuations as the target is approached. This difference 
arises because the transmission of forces within the gran- 
ular jet occurs via force chains. A small positional rear- 
rangement of the grains can cause a large rearrangement 
of how stresses distribute through the jet. Aside from 
these fluctuations, we nevertheless see good agreement. 
Granular jet impact, at leading order, is controlled by 
the same pressure profile as that governing the impact of 
a perfect fluid jet. 

The finding that dense granular impact onto a friction- 
less target is controlled, at leading order, by a particu- 
larly simple continuum behavior raises an open question: 
what are the necessary ingredients for impact to produce 
highly collimated ejecta? 

Another example of an unexpected emergence of a 
liquid-like response in a system of particles after impact 
comes from the relativistic heavy ion collider (RHIC). 
RHIC was designed to produce quark gluon plasma, a 
state of matter at a very high energy density obtainable 
only by means of heavy nucleon collisions. Collisions at 
RHIC did not produce an isotropic scattering pattern, 
as is found in collider experiments using lighter hadronic 
particles. Instead the ejecta is collimated as evidenced by 
elliptical momentum distribution in the ejecta. This col- 
limation has been interpreted as a signature that a liquid 
quark-gluon phase was created by the impact [I5j [16] . 
However, the presence of local equilibrium throughout 
the entire beam during impact at RHIC is currently not 
fully justified [17 j, and therefore coherent ejecta at RHIC 
might not be due to the material phase on impact. 

In dense granular jet impact we also find difficulty at- 
tributing collimation to a liquid phase, since the jet is 
far from equilibrium. This finding together with the fact 
that the presence or absence of a dead zone in the jet 
interior introduces only minor changes in the ejecta dy- 
namics suggests to us that impact generically induces a 
liquid-like ejecta response when the initial particle den- 
sity is sufficiently large. Collimation occurs regardless of 
what happens to the material within the high-pressure 
region. We have provided supporting evidence for this 
view by identifying the continuum limit controlling the 
ejecta formation as perfect fluid flow. In other words, 
the observed liquid-like response can be quantitatively 
accounted for by considering the mechanics of continuum 
impact alone. It is not necessary to invoke the existence 
of a liquid phase. 

Our results also shed light on important technologi- 
cal processes, such as blast cleaning, in which a stream 
of glass beads are accelerated towards engine blades in 
order to clean them. More broadly, the discovery that 
impact onto a frictionless target simplifies dense granu- 



lar flow into that of a perfect fluid allows us probe more 
realistic impact conditions analytically, via perturbative 
expansions, thereby potentially providing a more rigor- 
ous understanding of how dense granular flow organizes 
during impact. 

In conclusion, we demonstrate here that the impact of 
a dense granular jet onto a frictionless target is controlled 
by Euler flow. Predictions from an exact solution for 
perfect fluid impact quantitatively reproduce the veloc- 
ity and pressure fields found numerically. Discrepancies 
between the two solutions can be understood largely as a 
consequence of the fact that the perfect fluid impact so- 
lution conserves energy while dissipative processes, such 
as inelasticity and friction, are present within the gran- 
ular jet. This correspondence is important because it 
demonstrates that a liquid-like dynamical response, the 
ejection of a collimated sheet after impact, owes its exis- 
tence to incompressibility and momentum conservation, 
i.e. Eq. 0. The existence of a liquid phase is unneces- 
sary. 
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